Photoproduction of ix + n meson pairs on the proton 
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The exclusive reaction jp — > pn + -K~ was studied in the photon energy range 3.0 - 3.8 GeV and 
momentum transfer range 0.4 < -t< 1.0 GeV 2 . Data were collected with the CLAS detector at 
the Thomas Jefferson National Accelerator Facility. In this kinematic range the integrated lumi- 
nosity was about 20 pb _1 . The reaction was isolated by detecting the tv + and proton in CLAS, 
and reconstructing the n~ via the missing-mass technique. Moments of the di-pion decay angu- 
lar distributions were derived from the experimental data. Differential cross sections for the S, 
P, and D-waves in the M 1T + 7T - mass range 0.4 — 1.4 GeV were derived performing a partial wave 
expansion of the extracted moments. Besides the dominant contribution of the p(770) meson in 
the P-wave, evidence for the /o(980) and the /2(1270) mesons was found in the S and D-waves, 
respectively. The differential production cross sections da/dt for individual waves in the mass range 
of the above-mentioned mesons were extracted. This is the first time the /o(980) has been measured 
in a photoproduction experiment. 



I. INTRODUCTION 



The two pion channel offers the possibility of investi- 
gating various aspects of the meson resonance spectrum. 
It couples to the scalar-isoscalar channel that contains 
the tr, /o(980) and possibly a few more resonances with 
masses below 2 GeV. It is the main decay mode of the 
lowest isoscalar-tensor /2(1270) resonance and it is the 
only decay mode of the isovector-vector resonance, the 
p(770). Among all these, the p-meson is by far the most 
prominent and most extensively studied, both from the 
point of view of its production mechanisms and its inter- 
nal properties. Nowadays the other resonances too are 
subjects of extensive theoretical and experimental inves- 
tigation. The a meson is now established with pole mass 
and width determined with good accuracy [1-3]. How- 
ever, its microscopic structure seems to be quite different 
from that of the p and it is the subject of theoretical de- 
bate [4]. The /o(980) is even a more enigmatic state: its 
experimental determination is complicated by its proxim- 
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ity to the K K threshold, and its QCD nature still awaits 
an explanation [5]. Finally, the /2(1270) has been repre- 
sented so far as a Breit-Wigner resonance [2] and appears 
to fit well into the quark model spectrum [6] . 

In this paper we focus on the scalar sector, using the 
p meson as a benchmark for the analysis procedure. The 
KK channel from the same data set is currently being 
analyzed and in the near future a coupled-channels anal- 
ysis will provide further constraints on the extraction of 
the meson properties. 

For a long time most of our knowledge on the scalar 
meson spectrum was obtained from hadron-induced re- 
actions, 77 collisions and studying the decays of vari- 
ous mesons, e.g. </>, J/*, D and B. Very few studies 
were attempted with electromagnetic probes, in particu- 
lar real photons, since their production cross sections are 
relatively small compared to the dominant production 
of vector mesons. On one hand, through vector meson 
dominance, the photon can be effectively described as a 
virtual vector meson. On the other hand, quark-hadron 
duality and the point-like-nature of the photon coupling 
make it possible to describe photo-hadron interactions at 
the QCD level. Recently, high-intensity and high-quality 
tagged-photon beams, as the one available at JLab, have 
opened a new window into this field. 

In photoproduction processes, information about the 
S'-wave strength can be extracted by performing a par- 
tial wave analysis. Angular distributions of photopro- 
duced mesons and related obscrvablcs, such as the mo- 
ments of the angular distributions and the density matrix 
elements, are the most effective tools to look for inter- 
ference patterns. An interference between the .S'-wave 
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and the dominant P-wave was discovered in the moment 
analysis of K + K~ photoproduction on hydrogen, ana- 
lyzing the data collected in the experiments performed 
at DESY [7] and Daresbury [8]. In two-pion production 
experiments, such as reported in Rcfs. [9-11], moments 
and density matrix elements were used to analyze the 
properties of helicity amplitudes describing the photo- 
production process. Unfortunately, only the dominant 
spin-1 partial wave of the tt + tt~ pair was taken into ac- 
count. No attempt to obtain information about the S- 
wave amplitude was made. More recently, the HERMES 
experiment at DESY [12] investigated the interference of 
the P-wave in the w + ir~ system with the S and D- waves 
in the n + ir~ electroproduction process, and showed that 
such interference effects are measurable. The large pho- 
ton virtuality Q 2 >3 GeV 2 is, however, a crucial factor 
that distinguishes this analysis from the photoproduction 
analysis [9, 10]. 

Theoretical models for ir + ir~ photoproduction have 
been investigated in a series of articles. A very success- 
ful approach is the one by Soding [13] and its numerous 
modifications [14-17]. These models were able to de- 
scribe the shift of the maximum of the ir + n~ effective 
mass distribution with respect to the nominal p mass 
and the asymmetric shape observed in SLAC [9, 10] and 
DESY [11, 18] data. These properties are attributed to 
the interference of the dominating diffractivc p meson 
production, with its subsequent decay into 7r + 7r~, with 
the amplitudes corresponding to Drell-type diagrams in 
which the photon dissociates into n + and 7r _ , and one 
of the pions is elastically scattered off the proton. More 
recently, Gomez Tejedor and Oset [19] applied an effec- 
tive Lagrangian to construct the photoproduction am- 
plitudes. Their approach is limited to photon energies 
below 800 McV and effective masses M wn smaller than 
1 GeV. A two-stage approach for the tt + it~ 5-wave pho- 
toproduction was proposed in the model of Ref. [20]. 
First, a set of Born amplitudes, corresponding to pho- 
toproduction of 7r + 7r~, 7r°7r°, K + K~ and K°K° pairs 
is calculated. Then the photoproduccd meson pairs are 
subject to final-state interactions resulting in the 7r + 7r~ 
system [21-24]. The coupled-channels calculations were 
separately performed for all isospin / components of the 
transition matrix. Thus the 5-wave amplitudes in that 
model account for the existence of the isoscalar a, /o(980) 
and /o(1500), and the isovector ao(980) and ao(1450) res- 
onances. The coupling of the KK isovector channel with 
the nr] amplitude is described in Ref. [25]. 

All theoretical approaches described above do not con- 
sider explicitly the s-channel production of baryon res- 
onances contributing to the jmir final state. Data from 
Refs. [9, 11, 18], as well as from more recent experimen- 
tal studies [26] , indicate that the contribution of baryon 
resonances, such as A ++ and A , dominate at lower inci- 
dent photon energies (below 2 GeV). Furthermore, data 
obtained with the SAPHIR detector at ELSA for pho- 
ton energies between 0.5 GeV and 2.6 GeV show that 
the contribution of baryonic resonances to the 7r + p and 



7r~p mass distributions gradually decreases with photon 
energy. 

In this paper we review the results of the analysis 
of 7r + 7r~ photoproduction in the photon energy range 
3.0 - 3.8 GeV and momentum transfer squared —t be- 
tween 0.4 GeV 2 and 1 GeV 2 , where the di-pion effec- 
tive mass Mn-Tr varies from 0.4 GeV to 1.4 GeV. The 
main results were previously reported in Ref. [27]. We 
are not aware of any previous evidence of scalar mesons, 
in particular of the /o(980), in photoproduction of pion 
pairs. This effective mass region is dominated by the 
production of the p(770) resonance in the P-wave. From 
other experiments, such as pion-nuclcon collisions -K~p — > 
ir + ir~n [28, 29] or nucleon-antinucleon annihilation [30], 
there is some evidence that resonant states are formed 
in the 5-wave. These resonances have been neglected in 
previous experimental analyses of tt + it~ photoproduc- 
tion and, to our knowledge, the current analysis is the 
first one that explicitly takes into account the possibility 
that the S-wave is produced in the tt + tt~ system. 

In the following, some details are given on the exper- 
iment and data analysis (Sec. II), on the extraction of 
the angular moments of the di-pion system (Sec. Ill), 
and the fit of the moments using a dispersion relation 
(Sec. IV). Results of the partial wave analysis (differ- 
ential cross section for each partial wave and the spin 
density matrix elements) and the physics interpretation 
are reported in Sec. V. 



II. EXPERIMENTAL PROCEDURES AND 
DATA ANALYSIS 

A. The photon beam and the target 

The measurement was performed using the CLAS 
detector [31] in Hall B at Jefferson Lab with a 
bremsstrahlung photon beam produced by a continuous 
60-nA electron beam of energy Eq = 4.02 GeV imping- 
ing on a gold foil of thickness 8 x 10 -5 radiation lengths. 
A bremsstrahlung tagging system [32] with a photon en- 
ergy resolution of 0.1% Eq was used to tag photons in 
the energy range from 1.6 GeV to a maximum energy 
of 3.8 GeV. In this analysis only the high-energy part 
of the photon spectrum, ranging from 3.0 to 3.8 GeV, 
was used. e + e~ pairs produced by the interaction of 
the photon beam on a thin gold foil were used to contin- 
uously monitor the photon flux during the experiment. 
Absolute normalization was obtained by comparing the 
e + e~ pair rate with the photon flux measured by a total 
absorption lead-glass counter in dedicated low-intensity 
runs. The energy calibration of the Hall-B tagger system 
was performed both by a direct measurement of the e + e~ 
pairs produced by the incoming photons [33] and by ap- 
plying an over-constrained kinematic fit to the reaction 
IP ~ * ptt + it~ , where all particles in the final state were 
detected in CLAS [34]. The quality of the calibrations 
was checked by looking at the mass of known particles, 
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as well as their dependence on other kinematic variables 
(photon energy, detected particle momenta and angles). 

The target cell, a Mylar cylinder 4 cm in diameter 
and 40-cm long, was filled by liquid hydrogen at 20.4 K. 
The luminosity was obtained as the product of the tar- 
get density, target length and the incoming photon flux 
corrected for data-acquisition dead time. The overall sys- 
tematic uncertainty on the run luminosity was estimated 
to be in the range of 10%, dominated by the uncertainties 
on the photon flux. 



less for low-energy negative hadrons, which can be lost 
at forward angles because their paths are bent toward 
the beam line and out of the acceptance by the toroidal 
field. Coincidences between the photon tagger and the 
CLAS detector triggered the recording of the events. The 
trigger in CLAS required a coincidence between the TOF 
and the ST in at least two sectors, in order to select 
reactions with at least two charged particles in the final 
state. An integrated luminosity of 70 pb _1 (~ 20 pb _1 
in the range 3.0< E 1 <3.8 GeV) was accumulated in 50 
days of running in 2004. 




M x i (yp^p7i + X) (GeV ) 

FIG. 1: Missing mass squared for the reaction jp — > pir + X 
and the ir~ peak. The shaded area indicates the retained 
events. 
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FIG. 2: Two dimensional plot of the invariant masses ob- 
tained combining pairs of particles of the exclusive reaction 
7P — > pn + ir~ . 



B. The CLAS detector 

Outgoing hadrons were detected in the CLAS spec- 
trometer. Momentum information for charged particles 
was obtained via tracking through three regions of multi- 
wire drift chambers [35] within a toroidal magnetic field 
(~ 0.5 T) generated by six superconducting coils. The 
polarity of the field was set to bend the positive particles 
away from the beam line into the acceptance of the de- 
tector. Time-of-flight scintillators (TOF) were used for 
charged hadron identification [36]. The interaction time 
between the incoming photon and the target was mea- 
sured by the start counter (ST) [37]. This is made of 24 
strips of 2.2-mm thick plastic scintillator surrounding the 
hydrogen cell with a single-ended PMT-based read-out. 
A time resolution of ~300 ps was achieved. 

The CLAS momentum resolution, cr p /p, ranges from 
0.5 to 1%, depending on the kinematics. The detector 
geometrical acceptance for each positive particle in the 
relevant kinematic region is about 40%. It is somewhat 



C. Data analysis and reaction identification 

The raw data were passed through the standard CLAS 
reconstruction software to determine the four-momenta 
of detected particles. In this phase of the analysis, cor- 
rections were applied to account for the energy loss of 
charged particles in the target and surrounding materi- 
als, misalignments of the drift chamber's positions, and 
uncertainties in the value of the toroidal magnetic field. 

The reaction jp — > pi: + Tr~ was isolated detecting the 
proton and the ir + in the CLAS spectrometer, while the 
7r~ was reconstructed from the four-momenta of the de- 
tected particles by using the missing-mass technique. In 
this way the exclusivity of the reaction is ensured, keep- 
ing the contamination from the multi-pion background 
to a minimum. Figure 1 shows the tt~ missing mass 
squared. The background below the missing pion peak 
appears as a smooth contribution in the 7T7t invariant 
mass without creating narrow structures. 
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FIG. 3: Invariant masses obtained combining pairs of particles 
of the exclusive reaction yp — > pir + n~ . Upper panel M^^,-; 
lower panel left M p7r +; lower panel right M p7T - . Spectra are 
not corrected for the detector acceptance. 



To avoid edge regions in the detector acceptance, only 
events within a fiducial volume were retained in this anal- 
ysis. In the laboratory reference system, cuts were de- 
fined for the minimum hadron momentum (p pro ton > 
0.32 GeV and p^+ > 0.125 GeV), and the minimum 
and maximum azimuthal angles (0 proton ^+ > 10° and 
6 V + < 120°). The fiducial cuts were defined comparing 
in detail the experimental data distributions with the re- 
sults of the detector simulation. The minimum momen- 
tum cuts were tuned for different hadrons to take into 
account the energy loss by ionization of the particles. 



After all cuts, 41M events were identified as produced 
in the exclusive reaction yp — > p7r + 7r~. The other event 
topologies, with at least two hadrons in the final state 
(p7r _ , 7r + 7r~, p7r + 7r~), were not used since in the kine- 
matics of interest for this analysis (— t < 1 GeV 2 ), the 
collected data are about one order of magnitude less due 
to the detector acceptance. Figures 2 and 3 show the 
invariant mass spectra of the different combinations of 
particles in the final state. The p(770) dominates the 7T7t 
spectrum and the A(1232) ++ peak is clearly visible in 
the pir + invariant mass. Figure 2 shows a small overlap 
between the A(1232) ++ and the irir spectrum. Bary- 
onic resonances in the p7r~ invariant mass spectrum are 
less pronounced. It has to be noted that the projection 
of the baryon resonance peaks in the 7T7t spectrum re- 
sults in a smooth contribution and cannot create narrow 
structures. The effect of this background was extensively 
studied as discussed in Sec.VC. 



III. MOMENTS OF THE DI-PION ANGULAR 
DISTRIBUTION 

In this section we consider the analysis of moments of 
the di-pion angular distribution defined as: 



(Y LM )(Ey,t,M„„) 



47t / dfL 



dtdMmrdQn 



(1) 

where da is the differential cross section (in momentum 
transfer t and di-pion invariant mass M W7r ), Yj,M are 
spherical harmonic functions of degree L and order M, 
and fijr = (On, </>„•) are the polar and azimuthal angles of 
the 7r + flight direction in the ir + ir~ helicity rest frame. 
For the definition of the angles in the di-pion system we 
follow the convention of Ref. [9]. It follows from Eq. 1 
that, for a given Ey, t and di-pion mass M ww , (Voo) corre- 
sponds to the di-pion production differential cross section 
da/dtdM^. 

There are many advantages in defining and analyz- 
ing moments rather than proceeding via a direct partial 
wave fit of the angular distributions. Moments can be ex- 
pressed as bi-linear in terms of the partial waves and, de- 
pending on the particular combination of L and M, show 
specific sensitivity to a particular subset of them. In ad- 
dition, they can be directly and unambiguously derived 
from the data, allowing for a quantitative comparison 
to the same observables calculated in specific theoretical 
models. 

Extraction of moments requires that the measured an- 
gular distribution is corrected by the detector acceptance. 
We studied three methods for implementing acceptance 
corrections. In the first two methods, the moments were 
expanded in a model-independent way in a set of basis 
functions and, after weighting with Monte Carlo events, 
they were compared to the data by maximizing a likeli- 
hood function. The first of these two parametrizes the 
theory in terms of simplified amplitudes, while the sec- 
ond uses directly moments as defined above. The ap- 
proximations in these methods have to do with the choice 
of the basis and depend on the number of basis functions 
used. The systematic effect of such truncations was stud- 
ied and the main results are reported below. In the last 
method, data and Monte Carlo were binned in all kinc- 
matical variables. The data were then corrected by the 
acceptance defined as the ratio of reconstructed over gen- 
erated Monte Carlo events in that bin. Since it was found 
to be not reliable in bins where the acceptance was small 
or vanishing, this method was only used as a check of the 
others and was not included in the final determination of 
the experimental moments. 



A. Detector efficiency 

The CLAS detection efficiency for the reaction jp — > 
p7r + 7r~ was obtained by means of detailed Monte Carlo 
studies, which included knowledge of the full detector 
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geometry and a realistic response to traversing particles. 
Events were generated according to three-particle phase 
space with a brcmsstrahlung photon energy spectrum. 
A total of 4 billion events were generated in the energy 
range 3.0 GeV < E 1 < 3.8 GeV and covered the allowed 
kinematic range in —t and M ww . About 700M events 
were reconstructed in the M Tlr and —t ranges of inter- 
est (0.4 GeV < M„ n < 1.4 GeV , 0.1 GeV 2 < -t < 1.0 
GeV 2 ). This corresponds to more than fifteen times the 
statistics collected in the experiment, thereby introduc- 
ing a negligible statistical uncertainty with respect to the 
statistical uncertainty of the data. 



B. Extraction of the moments via likelihood fit of 
experimental data 

Moments were derived from the data using detec- 
tor efficiency-corrected fitting functions. As mentioned 
above, the expected theoretical yield was parametrized in 
terms of appropriate physics functions: production am- 
plitudes in one case and moments of the cross section in 
the other. The theoretical expectation, after correction 
for acceptance, was compared to the experimental yield. 
Parameters were extracted by maximizing a likelihood 
function defined as: 



V{Ta)I(T a ) 



JdT V (T)I(T) 



(2) 



Here a represents a data event, n = AN is the number 
of data events in a given (-Ely, —t, M„) bin (i.e. the fit is 
done independently in each bin) , r a represents the set of 
kinematical variables of the a th event, r?(r a ) is the corre- 
sponding acceptance derived by Monte Carlo simulations 
and I(r a ) is the theoretical function representing the ex- 
pected event distribution. The measure dr includes the 
phase space factor and the likelihood function is normal- 
ized to the expected number of events in the bin 



n= dTT](r)I(T). 



(3) 



The advantage of this approach lies in avoiding bin- 
ning the data and the large uncertainties related to the 
corrections in regions of CLAS with vanishing efficien- 
cies. Comparison of the results of the two different 
parametrizations allows one to estimate the systematic 
uncertainty related to the procedure. In the following, 
we describe the two approaches in more detail. 



1. Parametrization with amplitudes 
The expected theoretical yield in each bin is described 



I(r a ) = 4tt 



^ ^ a LM (E 11 ~t 1 M^)Y LM (Q. 7r ] 



L=0 M=-L 



•(4) 



This parametrization has the benefit that the intensity 
function I(r a ) is by construction positive. However, it 
can lead to ambiguous results since it has more parame- 
ters than can be determined from the data. In addition, 
for practical reasons, the parametrization involves a cut- 
off, L max , in the maximum number of amplitudes. For a 
specific choice of L max , the number of fit parameters is 
given by 2(L max + l) 2 . We also note that these ampli- 
tudes are not the same as the partial wave amplitudes in 
the usual sense of a di-pion photoproduction amplitude, 
since the latter depends on the nucleon and photon spins. 

After removing the irrelevant constants, the fit is per- 
formed minimizing the function: 

AN 

In £ = -J2^V(Ta)I(T a ) (5) 

a=l 

+ AVln ^ a* L , M ,a LM ^L'M';L,M, 

L'M'-LM 

where we have introduced the rescaled amplitudes 
aLM^Ey, —t, M^n) defined by: 



(6) 



and the acceptance matrix ty(E~, — t, M nv ) was com- 
puted using Monte Carlo events as: 



4tt 



AN R<!C 



AN 



Gen 



J2 Y£, m ,(SK)Ylm(SK), (7) 



where ANc en and ANn ec are the number of generated 
and reconstructed events, respectively. 

Fits were done using MINUIT with the analytical ex- 
pression for the gradient, and using the SIMPLEX pro- 
cedure followed by MIGRAD [38]. After each fit, the 
covariance matrix was checked and if it was not positive 
definite, the fit was restarted with random input param- 
eters. At the end, the uncertainties were computed from 
the full covariance matrix. 



2. Parametrization with moments 

The expected theoretical yield in each (E 1: — t, M W7r ) 
bin is described as: 



E 

L=0 M=0 



(8) 



The parametrization in terms of the moments directly 
gives the quantities we are interested in (moments 
(Ylm))- However, the fit has to be restricted to make 
sure the intensity is positive. As in the amplitude 
parametrization, a cutoff L max in the maximum num- 
ber of moments has to be used. The number of fit pa- 
rameters is given by (L max + l){L max + 2)/2. As L max 
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increases, moments with L close to L max show a signifi- 
cant variation, while moments with the lowest L remain 
unchanged. 

The expected (acceptance-corrected) distribution is 
then given by: 



Z(r„) = V^J2 ^ LM ReF iA r(Or)] (Y 



LM 



(9) 



L.M 



The function to be minimized with respect to (Ylm) (L > 
0) is then given by: 



AW 



-21n£ = -2^ In J(r a ), 



(10) 



with the coefficients rjLAiiE^, — t, M n7V ) computed using 
Monte Carlo events 



-in 



AN, 



Gen 



AN Rec 

E 



ReY LM (Sli) 



(11) 



where — 1 for L = and 1/2 for all other (LM). 
For L rnax < 4, the results are similar to what was ob- 
tained with the previous method, showing the same sta- 
bility against L max truncation and a similar goodness of 
the fit. To check the sensitivity of the likelihood fit to 
the parameter initialization, moments were extracted in 
three different ways: 1) using random initialization for all 
parameters; 2) fixing the parameters up to L = 2 to the 
ones obtained from a fit with L max = 2, and randomly 
initializing the others; 3) starting with parameters ob- 
tained in 2) and then releasing all parameters. The three 
different methods gave consistent results and the differ- 
ence of moments obtained using the different procedures 
was used to evaluate the systematic uncertainty related 
to the fit procedure. 



> 



V 



> 
1 



V 



> 
o 



V 




-2 



-4 



1.2 1.4 

M w (GeV) 





i 


»" 



0.6 



0.8 



1.2 1.4 
MTGeV) 




1.2 1.4 
M^GeV) 

FIG. 4: Moments of the di-pion angular distribution in 3.2 < 
£ 7 < 3.4 GeV and -t = 0.45±0.05 GeV 2 (black), -t = 0.65± 



0.05 GeV (red) and 



0.95 ± 0.05 GeV (blue). Error 



bars include both statistical and systematic uncertainties as 
explained in the text. 



3. Methods comparison and final results 

Moments derived by the different procedures agreed 
qualitatively. The most stable results were obtained by 
using the first parametrization, although we do find oc- 
casionally large bin-to-bin fluctuations. However, there 
are no a priori reasons to prefer one of the two methods 
and we consider the discrepancies between the fit results 
as a good estimate of the systematic uncertainty associ- 
ated with the moments extraction. The final results are 
given as the average of the first method (parametrization 
with amplitudes) and the second method (parametriza- 
tion with moments) with the three fit initializations: 



Y final 



E 

i—l, 4 Methods 



(12) 



where Y stands for (Ylm){E 1 , t, M OT7r ). 

The total uncertainty on the final moments was eval- 
uated adding in quadrature the statistical uncertainty, 



SYminuit as given by MINUIT, and two systematic un- 
certainty contributions: SY sys t fa related to the moment 
extraction procedure, and SY syst norm, the systematic un- 
certainty associated with the photon flux normalization 
(see Sec. II). 



SY, 



final 



— \l &Y MINUIT + SY syst j it + 5Yg yst norm (13) 



with: 



8Y syst fa 

i=l,4 Methods 
SY S yst norm 10% ' ^ 'final 



\ 



(14) 



(15) 



For most of the data points, the systematic uncertain- 
ties dominate over the statistical uncertainty. Samples of 
the final experimental moments are shown in Figs. 4, 5, 6, 
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FIG. 



5: Moments of the di-pion angular distribution in 3.2 < 

2 



E 1 < 3.4 GeV and -t = 0.45±0.05 GeV 2 (black), -t = 0.65± 
0.05 GeV 2 (red) and -t = 0.95 ± 0.05 GeV 2 (blue). Error 
bars include both statistical and systematic uncertainties as 
explained in the text. 



FIG. 6: Moments of the di-pion angular distribution in 3.2 < 
£ 7 < 3.4 GeV and -t = 0.45±0.05 GeV 2 (black), -t = 0.65± 
0.05 GeV 2 (red) and -t = 0.95 ± 0.05 GeV 2 (blue). Error 
bars include both statistical and systematic uncertainties as 
explained in the text. 



and 7. The whole set of moments resulting from this 
analysis is available at the Jefferson Lab [39] and the 
Durham [40] databases. 

As a check of the whole procedure, the differential cross 
section da/dt for the 77? — > pp(770) meson has been ex- 
tracted by fitting the (loo) moment in each —t bin with a 
Breit-Wigner plus a first-order polynomial background. 
The agreement within the quoted uncertainties with a 
previous CLAS measurement [41], as well as the world 
data [11], gives us confidence in the analysis procedure. 

IV. PARTIAL WAVE ANALYSIS 

In the previous section we discussed how moments of 
the angular distribution of the 7r + 7r _ system, (Ylm), 
were extracted from the data in each bin in photon en- 
ergy, momentum transfer and di-pion mass. In this sec- 
tion we describe how partial waves were parametrized 



and extracted by fitting the experimental moments. 

Moments can be expressed as bi-linear in terms of the 
amplitudes aim = ai m (\, A', A 7 , E~ h t, M^^) with angular 
momentum I and z-projection m (in the chosen reference 
system m coincides with the helicity of the di-pion sys- 
tem) as: 

(Y L m) = ^2 C(l'm',lm,LM) x ai m aj, m ,, (16) 

I'm' ,/m, A, A' 

where C are Wigner's 3jm coefficients, A 7 is the helicity of 
the photon, and A and A' arc the initial and final nucleon 
helicity, respectively The explicit forms of the moments 
with L < 4 in terms of amplitudes with / = (S'-wave), 
/ = 1 (P-wave), I = 2 (£>-wave), and I = 3 (P-wave) are 
given in Appendix A. 
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FIG. 7: Moments of the di-pion angular distribution in 3.2 < 
E 1 < 3.4 GeV and -t = 0.45±0.05 GeV 2 (black), -t = 0.65± 
0.05 GeV 2 (red) and -t = 0.95 ± 0.05 GeV 2 (blue). Error 
bars include both statistical and systematic uncertainties as 
explained in the text. 



A. Helicities, isospin and coupled-channels 
dependence 

The photon helicity was restricted to A 7 = +1 since 
the other amplitudes are related by parity conservation. 
In addition, some approximations in the parametrization 
of the partial waves were adopted to reduce the number 
of free parameters in the fit and are discussed below. 

• The number of waves was reduced restricting the 
analysis to |m| < 1 since to = 2 is only possible 
for I > 2 (D and F waves), which are expected to 
be small in the mass range considered [9, 10]. In 
the chosen reference system, m coincides with the 
helicity of the di-pion system and, since we used as 
a reference the wave with A 7 = +1, the three val- 
ues of m have a simple interpretation in terms of 
helicity transfer from the photon to the 7T7r-system: 
m = +1 corresponds to the non-helicity-flip am- 



plitude (s-channcl helicity conserving) that is ex- 
pected to be dominant [10], while m = 0, —1 corre- 
spond to one and two units of helicity flip, respec- 
tively. In the case of the S'-wavc (I = m = 0), only 
one amplitude is considered. 

• The dependence on the nucleon helicity was sim- 
plified as follows. For a given l,m,E^,t set, there 
are four independent partial wave amplitudes cor- 
responding to the four combinations of initial and 
final nucleon helicity, A and A'. It is expected that 
the dominant amplitudes require no nucleon helic- 
ity flip [10]. Without nucleon polarization informa- 
tion it is not possible to extract all four amplitudes. 
Thus our strategy is to consider in the analysis only 
the dominant ones or to exploit possible relations 
among them. For example, in the Regge p and uj 
exchange model, the following relations are satisfied 

by the 5-wave amplitudes: (AA') = (++) = ( ) 

and (H — ) = —(—+), where ± corresponds to helic- 
ity ±1/2. More generally, by examining the exper- 
imental moments, we observe that the interference 
between the dominant P-wave, seen in the (Y21) 
moment in the p region, indicates that the P m =+i 
and the P m =o amplitudes are out of phase. For a 
single nucleon-helicity amplitude, this would imply 
a difference between the (Yn) and (Y 10 ) moments, 
arising primarily from the interference between the 
5*-wave and the P m=+ \ and P m =o waves, respec- 
tively, in the p region where the S amplitude does 
not vary substantially. The data suggests, however, 
that both (Y11) and (Y10) peak near the position 
of the p. A possible explanation for the behavior 
of the data is the following: the dominant P m= +\ 
amplitude may originate from the helicity-non-flip 
diffractive process and the P m =o amplitude from a 
nucleon-hclicity-flip vector exchange, which is also 
expected to contribute to the S'-wave production. 
This would also explain why the (Yn) and (Y w ) 
moments have comparable magnitudes. To accom- 
modate such behavior, at least two nucleon-helicity 
amplitudes are required. In addition, since strong 
interactions conserve isospin, it is convenient to 
write the 7T7T amplitudes in the isospin basis. Each 
amplitude was then expressed as a linear combi- 
nation of 7T7T amplitudes of fixed isospin / (with 
/ = 0,1,2). 

• The coupling of the 7T7t system to other chan- 
nels was taken into account introducing a multi- 
dimension channel space: for a given isospin / in 
the partial wave I, the amplitudes depend also on 
an index a that runs over different di-meson sys- 
tems. For example, a = 1 corresponds to 7T7t, a = 2 
to KK, a = 3 to 7777, etc. In the subsequent analy- 
sis we will restrict the channel space to include the 
7T7r and KK channels, which are the only channels 
relevant in the energy range considered. 

According to these considerations, the moments 
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FIG. 8: Fit result (black line) of the final experimental mo- 
ments (in red) for 3.2 < _E 7 < 3.4 GeV and 0.5 < —t < 
0.6 GeV 2 . The systematic uncertainty and fit uncertainty are 
added in quadrature and are shown by the gray band. 



where PV represents the principal value of the integral 
and p corresponds to the phase space term. In this ex- 
pression, / and Simj are matrices in the multi-channel 
space (7T7T, KK) 1 as mentioned above. Ni m j and Di m j 
can be written in terms of the scattering matrix of tt-k 
scattering, chosen to reproduce the known phase shifts, 
inelasticities [49, 50], and the isoscalar (I = S, D), isovec- 
tor (I = P, F) and isotensor (I = S, D) amplitudes in the 
range 0.4 GeV < </s < 1.4 GeV. Finally, the amplitude 
Q-im.i represents our ignorance about the production pro- 
cess. 

As a function of s = M„ n , ai m ,i have cuts for s > 4m| 
(right-hand cut) and for s < m 2 (left-hand cut). The 
left-hand cut reflects the nature of particle exchanges de- 
termining the 7T7T photoproduction amplitude, while the 
right-hand cut accounts for the final-state interactions 
of the produced pions. In Eq. 18, these discontinuities 
are taken into account by the functions Ni m j and Di m j, 
while ai m j(s) does not have singularities for s > 4m„ 
and can be expanded in a Taylor series: 



air, 



[A + Bs + Cs 2 + ■■■] [k] 



(19) 



with A, B, . . . being matrices of numerical coefficients to 
be determined by the simultaneous fit of the angular mo- 
ments defined in Eq. 16 and [k] = k l a 8 a ^ used to take into 
account the threshold behavior in the l-th partial wave. 
All amplitudes but the scalar-isoscalar are saturated by 
the 7T7T state. For the scalar-isoscalar amplitude, the K K 
channel was also included. In addition, to reduce sensi- 
tivity to the large energy behavior of the (tttTjKK) am- 
plitudes, the real part of the integral was subtracted and 
replaced by a polynomial in s, whose coefficients were 
also fitted. 



RESULTS 



were fitted to a set of amplitudes given by: 

for each l,m, \m\ < 1, with i = 1,2 corresponding 
to the nucleon helicity non-flip and helicity-ffip of 
one unit, isospin / = 0, 1, 2 and channel a. 



B. Amplitude parametrization 

For each helicity state of the target A, recoil nucleon 
A', and tttt system m, in a given E 1 and t bin, the cor- 
responding helicity amplitude ai m (s = M%„), was ex- 
pressed using a dispersion relation [42-48] as follows: 

ai m ,i(s) = i[I + Si m j(s)]a lm j(s) (18) 

•* 8th S S 



A. Fit of the moments 

Using the parametrization of the partial waves de- 
scribed in the previous section, we fitted all moments 
(Ylm) with L < 4 and M < 2 using amplitudes with 
I < 3 (up to F- waves). In Figs. 8, 9, 10, and 11 we 
present a sample of the fit results for E 1 = 3.3 ±0.1 GeV 
and 0.5 < |t| < 0.6 GeV 2 . 

To properly take into account the statistical and sys- 
tematic uncertainty contributions to the experimental 
moments described in Sec. Ill, the four sets of moments 
resulting from the different fit procedures were individ- 
ually fitted and the results were averaged, obtaining the 
central value shown by the black line in the figures. The 
error band, shown as a gray area, was calculated fol- 
lowing the same procedure adopted for the experimental 
moments. The final uncertainty was computed as the 
sum in quadrature of the statistical uncertainty of the fit 
and the two systematic uncertainty contributions. The 
first is related to the moment extraction procedure and 
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FIG. 9: Fit result (black line) of the final experimental mo- 
ments (in red) for 3.2 < _E 7 < 3.4 GeV and 0.5 < -t < 
0.6 GeV 2 . The systematic uncertainty and fit uncertainty are 
added in quadrature and are shown by the gray band. 



FIG. 10: Fit result (black line) of the final experimental mo- 
ments (in red) for 3.2 < E 1 < 3.4 GeV and 0.5 < -t < 
0.6 GeV 2 . The systematic uncertainty and fit uncertainty are 
added in quadrature and are shown by the gray band. 



is evaluated as the variance of the four fit results. The 
latter is associated with the photon flux normalization 
and is estimated to be 10%. The central values and un- 
certainties for all the observables of interest discussed in 
the following sections were derived from the fit results 
with the same procedure. 

The moment (loo), corresponding to the differential 
production cross section da / dtdM , shows the dominant 
p(770) meson peak. In the (Yio) and (Yn) moments, the 
contribution of the S'-wave is maximum and enters via 
interference with the P-wave. In particular the struc- 
ture at ~ 0.77 GeV in (Yn) is due to the interfer- 
ence of the S- wave with the dominant, helicity-non-flip 
wave, P m=+ i. In the (Yin) moment the same structure 
is due to the interference with the P m =o wave, which 
corresponds to one unit of hclicity flip. A second dip 
near M„ = 1 GeV is clearly visible and corresponds to 
the production of a resonance that we interpret as the 
Jo (980). 



B. Partial wave amplitudes 

The square of the magnitude of the S-, P-, D- and 
F- waves resulting from the fit, summed over the nucleon 
spin projections, is given by: 

hm = 2_j \ a im,%{E^, t, M T7r )| 2 . 

i=l,2 

(20) 

When summed over the di-pion helicity, this can be writ- 
ten as: 

h = ^ \aim,i{E 7 , t, M n7r )\ 2 , 

(21) 

where the sum is limited to m = — 1, 0, 1 for I > and to 
m = for I = 0. 

The resulting partial wave cross sections are shown in 
Figs. 12, 13, 14, and 15, for a selected photon energy 
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FIG. 11: Fit result (black line) of the final experimental mo- 
ments (in red) for 3.2 < P 7 < 3.4 GeV and 0.5 < —t < 
0.6 GeV 2 . The systematic uncertainty and fit uncertainty are 
added in quadrature and are shown by the gray band. 




FIG. 12: S-wave cross section derived by the fit in the 
3.2 < E 1 < 3.4 GeV and 0.5 < -t < 0.6 GeV 2 bin. The 
systematic and the fit uncertainties are added in quadrature 
and are shown by the gray band. 




1.4 

(GeV) 



FIG. 13: P-wave cross section derived by the fit in the 
3.2 < P 7 < 3.4 GeV and 0.5 < -t < 0.6 GeV 2 bin. Bot- 
tom plots: the same amplitudes for the three possible values 
of An-vr (from left to right -1,0 and +1). The systematic and 
fit uncertainties are added in quadrature and are shown by 
the gray band. 



and —t bin. The whole set of partial wave amplitudes 
resulting from this analysis is available at the Jefferson 
Lab [39] and the Durham [40] databases. 

As expected, the dominant contribution from the p me- 
son is clearly visible in the P-wave, whose contribution 
is about one order of magnitude larger than the other 
waves. In particular the main contribution comes from 
Iim=i,+ii corresponding to a non-helicity flip (s-channcl 
hclicity conserving) transition. In the S'-wavc, a strong 
interference pattern shows up around M X7r = 980 MeV, 
which reveals contributions from the Jo (980) production. 
The contribution from the /2(1270) tensor meson is ap- 
parent in the D-wave, while no clear structures are seen 
in the P-wavc. 



C. Systematic studies 

The error bands plotted in Figs. 12, 13, 14, and 15 
include the systematic uncertainties related to the mo- 
ment extraction and the photon flux normalization as 
discussed in Sec. Ill B 3. In addition, for the S'-wave, 
where the fo (980) contribution is strongly affected by in- 
terference, detailed systematic studies using both Monte 
Carlo and data were performed. 

In order to test the approximation introduced by the 
truncation to L max =A in the moment extraction, we first 
verified the fit was able to reproduce the experimental 



13 





M ra (GeV) 





















1 

1 



0.4 



1.4 



M^GeV) M^CGeV) (GeV) 



FIG. 14: As Fig. 13 for L>-wave. 
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FIG. 15: As Fig. 13 for F-wa,ve. 



distributions in the kinematic range of interest. Figure 16 
shows the comparison between data and fit results for 
the decay angles in the hclicity system with M wv in the 
/o(980) mass region (M wv = 0.985±0.01 GeV). Figure 17 
shows the same comparison for the invariant mass M p7r + 
when three different regions of M vv (M nv = 0.475 ± 
0.01 GeV, M„ n = 0.775 ± 0.01 GeV, M nn = 1.295 ± 
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FIG. 16: Pion angles in the -k + tt helicity rest frame for M nn 
in the / (980) mass region (M nn = 0.985 ± 0.01 GeV). Ex- 
perimental data are plotted in black and fit results in red. 



0.01 GeV) were selected. The good agreement proves 
the accuracy of the approximation. 

As a second check, we applied the fit to pseudo-data 
obtained with a realistic event generator, processed with 
the CLAS GEANT-based simulation package and ana- 
lyzed with the same procedure used for the data. Since 
the event generator was tuned to previous two-pion pho- 
toproduction measurements, it does not include any ex- 
plicit limitation on the number of waves. The recon- 



structed moments showed that, with the chosen L r 



all 



fits were capable of reproducing the generated moments 
up to M V7r ~ 1.1 GeV. Finally, we derived a quantitative 
estimate of the truncation effect on the 5-wave squared 
amplitude as follows. The results of a L max = 8 fit of the 
moments was used as input for a new Monte Carlo event 
generator. After being processed in the same way as dis- 
cussed above, pseudo-data were fitted with L max = 4 and 
the iS-wave amplitude was extracted. The difference be- 
tween the generated and the reconstructed partial wave 
cross section was found to be of the order of 25% that, 
added in quadrature to the other systematic uncertain- 
ties, was included in the gray band of Fig. 12. 

We also demonstrated that no structures similar to the 
narrow interference pattern we are interpreting as the ev- 
idence of the /o(980) were created by distortions induced 
by the CLAS acceptance. This check was performed gen- 
erating events after removing the /o(980) contribution, 
and verifying that no spurious structures appeared in the 
spectra after the full GEANT simulation and reconstruc- 
tion. 
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FIG. 18: 5-, P- and D-wave cross sections in the 3.4 < 
E-, < 3.6 GeV and 0.5 < -t < 0.6 GeV 2 bin. The gray 
and red bands show the resuits of the standard fit and of 
the fit performed adding a second-order polynomial to the 
partial wave expansion of the moments to account for the 
baryon resonance contributions. The width of the bands rep- 
resents the fit uncertainties only. Fit results are shown for a 
specific parametrization of the moments (second method, see 
Sec. IIIB3). 



In addition, the effects of baryon resonance contribu- 
tions to the di-pion mass spectrum were studied perform- 
ing the fit of the moments with the inclusion of an inco- 
herent background. In fact, the background in the di- 
pion mass spectrum introduced by the reflection of the 
baryon resonances is expected to be smooth and struc- 
tureless, contributing to all waves. Therefore this was 
parametrized as a second-order polynomial in M WT that 
was summed to the parametrization of the moments in 
terms of partial waves used in the standard fits. From 
this study we concluded that the background contribu- 
tion is small, smooth and does not affect the quality of 
the fit. The comparison of the fit results with and with- 
out the inclusion of this additional background indicates 
that the P-wave and the ,5- wave in the /o(980) region are 
only slightly affected, as shown in Fig. 18. On the con- 
trary, the low mass S-wave, corresponding to the <r(600) 
region, and the £>-wavc, corresponding to the /2(1270) 
region, show a significant variation and, therefore, a more 
complete analysis should be performed to extract reliable 
information in these mass ranges. A similar conclusion 
was drawn by comparing the analysis results excluding 



the A(1232), the dominant baryon resonance contribu- 
tion for this final state, with the cut M(pTi + ) > 1.4 GeV. 
A negligible effect was found on the rapid motion around 
the narrow / (980) meson, while a larger variation was 
observed at higher values of the M(7T7t) mass. 

To verify the stability of the fit of moments in the 
region of the Jo (980), the whole analysis was repeated 
reducing the M vv bin size from 10 to 5 MeV. The results 
obtained in the two cases were found to be consistent. 

As a final check, the sensitivity to the specific choice 
of the number of terms used in the Taylor expansion of 
the amplitudes a L (see Eq. 19) was tested performing the 
partial wave analysis fits both with a second- and fourth- 
order polynomial. The effect was found to be negligible 
compared to the other systematic uncertainties. 



D. The spin density matrix elements 

From the production amplitudes derived by the fit, we 
calculated the spin density matrix elements [51] for the 
P-wave and the interference between the S- and P-waves. 
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FIG. 19: Spin density matrix elements for the P-wave in the 
3.0 < E 1 < 3.2 GeV and 0.4 < -t < 0.5 GeV 2 bin. The 
black dots are data points from Ref. [9], taken in a similar 
kinematic bin (P 7 ~ 2.8 GeV and 0.02 < -t < 0.4 GeV 2 ). 



Some selected results are shown in Figs. 19, 20 and 
21. Since these observables do not depend on the photon 
flux normalization, the error bands do not include the 
10% uncertainty mentioned above. The whole set of spin 
density matrix elements resulting from this analysis is 
available at the Jefferson Lab [39] and the Durham [40] 
databases. 

Comparisons of our measurements at 3.0 < E y < 
3.2 GeV and 0.4 < -t < 0.5 GeV 2 with existing data 
from Rcfs. [9, 10] in a similar kinematic domain {E 1 ~ 2.8 
GeV and 0.02 < -t < 0.4 GeV 2 ) arc shown in Fig. 19. 
As expected, the two matrix elements pro and pn agree 
very well since they have a weak dependence on —t, while 
Poo shows a similar behavior, but with different values as 
it is more sensitive to the momentum transfer. If one 
compares the larger —t bins we measured, the differences 
increase, showing that extrapolating our data to lower 
—t would probably give good agreement with previous 
measurements. 




(GeV) 

FIG. 20: Spin density matrix elements for the P-wave in the 
3.2 < £L < 3.4 GeV and 0.5 < -t < 0.6 GeV 2 bin. 
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FIG. 21: Spin density matrix elements for the interference 
between S- and P-waves in the 3.2 < E 1 < 3.4 GeV and 
0.5 < -t < 0.6 GeV 2 bin. 



E. Differential cross sections 



As shown in Fig. 21, around — 980 MeV an in- 
terference pattern clearly shows up in the S-P wave in- 
terference term, corresponding to the contribution from 
the /o(980) meson. 



The differential cross sections [da / dt]i_ wave for indi- 
vidual waves and mass resonance regions were obtained 
integrating the corresponding amplitudes. The cross sec- 
tions in the mass regions of the /o(980), p, and /2(1270) 
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mesons were obtained integrating the S-, P- and D-waves 
in the mass ranges 0.98 ± 0.04 GeV, 0.4-1.2 GeV, and 
1.275 ± 0.185 GeV, respectively. These are shown in 
Figs. 22, 23 and 24 in the photon energy range 3.0- 
3.8 GeV. As mentioned previously, the P-wave is com- 
pletely dominated by the p meson production, and there- 
fore the integrated cross section can be directly compared 
to the world's data for the jp — > pp reaction [11, 41]. It 
should be noticed that the previous cross sections were 
evaluated without performing a partial wave analysis but 
fitting the mass-dependent cross section with a relativis- 
tic Breit-Wigner plus a smooth polynomial function to 
separate the resonance from the background. The good 
agreement shown in Fig. 23 gives confidence in the partial 
wave analysis. As expected, the 5-wave photoproduction 
is suppressed compared to the P-wave by more than an 
order of magnitude, reflecting the different mechanisms 
that lead to scalar and vector meson photoproduction: 
in Rcggc theory the latter is dominated by Pomeron ex- 
change, while the former is dominated by the exchange of 
reggeons that become suppressed as the energy increases. 



VI. 



SUMMARY 



In summary, we have performed a partial wave anal- 
ysis of the reaction jp — > pir + Tr~ in the photon en- 
ergy range 3.0-3.8 GeV and momentum transfer range 
—t = 0.4 — 1.0 GeV 2 . Moments of the di-pion angular 
distribution, defined as bi-linear functions of partial wave 
amplitudes, were fitted to the experimental data with an 
unbinncd likelihood procedure. Different parametriza- 
tion bases were used and detailed systematic checks were 
performed to insure the reliability of the analysis pro- 
cedure. We extracted moments (Ylm) with L < 4 and 
M < 2 using amplitudes with I < 3 (up to F- waves). Us- 
ing a dispersion relation, unitarity constraint, and phase 
shifts and inelasticities of tttt scattering, the production 
amplitudes were expressed in a simplified form, where 
the unknown part was expanded in a Taylor series. The 
coefficients were fitted to the experimental moments to 
extract the S-, P-, D-, and P- waves in the M M range 
0.4-1.4 GeV. 




t (GeV ) 



FIG. 22: Differential cross section da/dt for the S- wave in 
the Mtttt range 0.98 ± 0.04 GeV and photon energy range 
E-, = 3.0 - 3.8 GeV. 
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FIG. 23: Differential cross section da/dt for the P-wave in 
the Mtttt range 0.4-1.2 GeV and photon energy range E 1 — 
3.0 - 3.8 GeV. 



good agreement with previous measurements, and for the 
first time, the S — P interference term was extracted. 



The moment (loo) is dominated by the p(770) meson 
contribution in the P-wave, while the moments (Yio) 
and (iii) show contributions of the S'-wave through 
interference with the P-wave. The clear structure at 
M W7r ~ 1 GeV seen in such experimental moments and 
in the 5-wave amplitude is evidence of a resonance con- 
tribution that we interpret as the /o(980). This is the 
first observation of the /o(980) scalar meson in photopro- 
duction. A contribution from the /2(1270) tensor meson 
was observed in the D-wave, while no resonant structures 
were seen in the P-wave. The cross sections of individual 
partial waves in the mass range of the p(770), /o(980), 
and /2(1270) were computed. Finally, the spin density 
matrix elements for the P-wave were evaluated, finding 
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APPENDIX A 



t (GeV) 



FIG. 24: Differential cross section da/dt for the D-wave in 
the Mtttt range 1.275 ± 0.185 GeV and photon energy range 
E~ = 3.0 - 3.8 GeV. 



The explicit expressions for the moments, defined in 
Eq. 1 in terms of partial waves, given Eq. 4, truncated to 
the L = 3 (P) wave are given by, 



(Yoo) 
(Yio) 



\S\ 2 + |P_| 2 + |P | 2 + \P+\ 2 + \D-\ 2 + \D \ 2 + |P+| 2 + |F_| 2 + |P | 2 + \F+\ 2 



+ 



SP* + P S 
(2A 



[P-D*_ + P*P_ + P+P*_ + D+P* 



(P D* + D Q P* 



(Yu) = 



(Y: 



2(1/ 



(Y2l) 
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[D-Fl + P_P* + D+Fl + F+D* + ) 

-PS* - SP* + P+S* + sp+) 
3 



216 
280 



(D F* + FqDq) 



20 
9" 
70 



P D*_ - D_P* + P D* + + D+P* 
P P* - F-D* + P P + + F + D* Q ) 



'P-D* + D P1 - P+P* - P P+) 

{D-Fq + F D*_ - P+P* - F D* + ) 



SD* +D Q S 
^27 



(2\P \ 2 -\P_\ 2 -\P+\ 2 + \F_\ 2 + \F+\ 



/ 5 /90 
, 5 (P F* + F P*) + (|P_| 2 + \D + \ 2 ) + ]/^\Do\ 2 



35 
\\Fo\ 



P_F* + p_P* + P+P* + P + P;) 



16, „ |2 



i (SD*+ + D+S* - SD*_ - DS*) + (P P; + P+P * - P-P* - P P*) 
^ (P_P* + F P* - P+Pq - F P* + ) + ^Jt (P P; + F+P* - PoFl - F_P *) 
{D D* + + P+P* - P P* - D-D*) + ( Fo p; + F+F* - F F* - F_F*) 



-(p_p; + p + p: 



— {-f+f* - p_pr 

30 v + " +> 




( y P-F* + P+P* + P+P* + F_P*J 



-D-D* - D+D*_ 



/ 1 8 / 1 08 

(Y30) = SF* + F S* + y — (-P_P*_ - D_P* - P+P; - P+P + ) + ^ — (P P * + P P *) 

+ sj— (d_f* + p_p* + p+p; + p+p;) + y — (p f* + f p*) 
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1 / 18 

- {sf; + f+s* - sfi - F-sr) + ^ — (p+d* q + d p; - p_d* - d pi] 



(Y32) = 
+ 

(no = 
+ 

(Y42) = 



[P D*. + D+P* - P D* - D_P*) - 
[DoF; + P+D* - D F1 - P_D*) 
(-P+D* - D-Pl - P_D* + - D+P*) 



(D+F* + P q D; - D_P* - F Q D*_ 



-D+P* - P_D; - D_p; - P+D* ) 



(-P+P+ - P+P+ - P-Fl - P_P!) 



(P P * + PoP *) 



(-|d h 



ID J 



I 2 + IP-I 2 ) 



' A (P + ^* + F P; - P_P* - PoP^) + ^ ^ ( PoF ; + F+Po * _ p p* _ p_p *j 
S (DoD; + D+D* - D-Dq - D D*_) + (P P; + F+F* - F F1 - F-F*] 

~5~ /To" 

'_ (- P+ f*_ - p_p; - p_p; - p+p:) + y- (-d_d; - d+d* ) 
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It follows from Eq. 1 that the (loo) moment is normal- 
ized by the differential cross section via, 
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